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Abstract 

We investigate the “family” relationship of a possible scalar nonet com¬ 
posed of the ao(980), the /o(980) and the a and k type states found in 
recent treatments of vrTT and vrA scattering. We work in the effective La- 
grangian framework, starting from terms which yield “ideal mixing” accord¬ 
ing to Okubo’s original formulation. It is noted that there is another solution 
corresponding to dual ideal mixing which agrees with Jaffe’s picture of scalars 
as qqqq states rather than qq states. At the Lagrangian level there is no dif¬ 
ference in the formulation of the two cases (other than the numerical values 
of the coefficients). In order to agree with experiment, additional mass and 
coupling terms which break ideal mixing are included. The resulting model 
turns out to be closer to dual ideal mixing than to conventional ideal mix¬ 
ing; the scalar mixing angle is roughly —17° in a convention where dual ideal 
mixing is 0°. 
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I. INTRODUCTION 


Recently there has been renewed discussion |l|- [19| about evidence for low energy broad 


scalar resonances in the tttt and tiK scattering channels. In the approach on which 
the present paper is based, a need was found for a tttt resonance (a) at 560 MeV and a ttK 
resonance (k) around 900 MeV. That approach, motivated by the 1/Nc approximation 
to QCD, involves suitably regularized (near the poles) tree level diagrams computed from a 
chiral Lagrangian and containing resonances within the energy range of interest. Attention 
is focussed on the real parts which satisfy crossing symmetry but may in general violate the 
unitarity bounds. Then the unknown parameters (properties of the scalars) are adjusted to 
satisfy the unitarity bounds (i.e. to agree with experiment). In this way an approximate 
amplitude satisfying both crossing symmetry and unitarity is obtained. 

Similar results for the scalars have been obtained in different models |]^ although 
there is not unanimous agreement. These are, after all, attempts to go beyond the energy 
region where chiral perturbation theory can provide a practical systematic framework. 

Now if one accepts a light a and k and notes the existence of the isovector scalar ao(980) 
as well as the /o(980) there are exactly enough candidates to £11 up a nonet of scalars, all 
lying below 1 GeV. Presumably these are not the “conventional” p-wave quark-antiquark 
scalars but something different. It would then be necessary (see for example the discussion 
on page 355 of [^) to have an additional nonet of “conventional” heavier scalars. 

Most mesons fit nicely into a pattern where they have quantum numbers of quark- 
antiquark (qq) bound states with various orbital angular momenta. Furthermore, their 
masses and decays are (roughly) explained according to a nonet scheme, first proposed by 
Okubo known as “ideal mixing”. It has been widely recognized that the low-lying 
scalars (at least the well observed ao(980) and /o(980)) do not appear to fit this usual 
pattern. Hence Jaffe [Q proposed an attractive scheme, in the context of the MIT bag 
model [^, in which the light scalars are taken to have a qqqq quark structure (and zero 
relative orbital angular momenta). Other models explaining light scalars as “meson-meson” 
molecules or as due to unitarity corrections related to strong meson-meson interactions 
also involve four quarks at the microscopic level and may possibly be related. 

Our concern in the present paper is to study the nonet structure of the light scalars 
based on the approach of H- H. There, an effective chiral Lagrangian treatment was used. 


In such a treatment, only the SU{3) flavor properties of the scalars are relevant |^. At 
this level, one would not expect any difference in the formulation of our model since both 
Okubo’s model and Jaffe’s model use nonets with the same SU{3) flavor transformation 
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properties. In fact, we shall show (in Section II) that the effective Lagrangian dehning ideal 
mixing in Okubo’s scheme has two “solutions”. The one he choses explains the light vector 
mesons with a natural quark-antiquark structure. The other solution is identical to Jaffe’s 
model of the scalars. We note that it may be formally regarded as having a dual-quark 
dual-antiquark structure, where the dual quark is actually an anti-diquark. 

The initial appearance is that the four masses of the light nonet candidates obey the 
ordering relation [Eq. (|2.9|) below] of the dual ideal mixing picture but not the more stringent 
requirement of this picture Eq. (|2.4|) . Furthermore the decay /o(980) —>■ tttt is experimentally 
observed but is predicted to vanish according to ideal mixing. Thus, it is necessary to 
consider some corrections to the ideal mixing model. When such correction terms are added 
[to yield a structure like Eq. ( ^.1(J| )] the new model actually displays two different solutions 
for the particle eigenstates corresponding to a given scalar mass spectrum (see the discussion 
in Section III) so it becomes unclear as to whether the ordinary or the dual ideal mixing 
picture is more nearly correct. In order to resolve this question the predictions for the 
scalar-pseudoscalar-pseudoscalar coupling constants are hrst computed for each of these two 
solutions. The hve coupling constants needed for ttK scattering are found to depend on 
only two parameters - A and B in Eq. (|3.8|) . Then (see Section IV) the ttK scattering 
is recalculated taking these two parameters as quantities to be £t. However it turns out 
that both solutions yield equally probable hts to the tiK scattering amplitudes. Finally, 
the question is resolved by noting that only one of the two solution sets gives results which 
could be compatible with the previous ^ tttt scattering analysis and with the /o(980) —>• tttt 
decay rate. 

The favored solution is characterized by a scalar cr — fo mixing angle which is closer to the 
dual form of ideal mixing than to the usual form. Using a convention [see Eq.( |3.6|) ] where 
an angle = 0 means dual ideal mixing and |0s| = | means conventional ideal mixing, the 
favored solution has Og ^ —17°. It should be noted that this result is based on an analysis 
of scalar coupling constants which are related to each other “kinematically” but which are 
related to experiment through “dynamical” models of ttK and tttt scattering. 

Some technical details are put in three Appendixes. Appendix A contains a brief discus¬ 
sion of some key features of the qqqq scalars as expected in the quark model. Appendix B 
shows how the needed terms of the Lagrangian including the scalar nonet may be presented in 
chiral covariant form. Finally Appendix C contains a list of the various scalar-pseudoscalar- 
pseudoscalar coupling constants and their relations to the parameters of our Lagrangian and 
to the scalar and pseudoscalar mixing angles. 
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II. SCALAR NONET MASSES 


For orientation, it may be useful to start off by paraphrasing Okubo’s classic discussion 

In our case the held will have = O’*" rather than = 1~ as in the original case. The 
notation is such that a lower index transforms under havor SU{3) in the same way as a quark 
while an upper index transforms in the same way as an antiquark. In this discussion it is not 
strictly necessary to mention the quark substructure of N - only its havor transformation 
property will be of relevance. This lack of specihcity turns out to be an advantage for our 
present purpose. 

The “ideal mixing” model may be dehned by the following mass terms of an ehective 
Lagrangian density: 


^ of the “ideal mixing” of a meson nonet held, which we denote as the 3x3 matrix N^{x). 


C^ass = -aTi{NN) - hTi{NNM), 


( 2 . 1 ) 


where a and b are real constants while M. is the “spurion matrix” M. = diag(l,l,a;) , x 
being the ratio of strange to non-strange quark masses in the usual interpretation. Iso-spin 
invariance is being assumed. The names of the scalar particles with non-trivial quantum 
numbers are: 


N 


Nl Oq k"*" 
Oo iV| 

K~ iV| 


( 2 . 2 ) 


with Oq = {Nl — A|)/a/ 2. There are two iso-singlet states: the combination {Nl + iV| -|- 
N^)/'/3 is an SU{3) singlet while {Nl + iV| — 2N^)/\/6 belongs to an SU{3) octet. These 
will in general mix with each other when SU{3) is broken. Diagonalizing the helds in Eq. 
0 yields the diagonal (ideally mixed) states {Nl + A|)/a/ 2 and A|. 

Now it is easy to read oh the particle masses from Eq. (0 in terms of a, b and x. This 
information is conveniently described by the two sum rules: 

9. ^ 9fNl+N^\ , , 

m^(ao)=m f—-^=—j, (2.3) 


(oo) — {k) = m? {k) — ■ (2-4) 

There are two characteristically diherent kinds of solutions, depending on whether both sides 
of Eq. ( |2.4|) are positive or negative. Okubo’s original scheme amounts to the choice that 
both sides of Eq. (|2.4|) are negative. Then 
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(2.5) 




(iV|) > {k) > (ao) = ) 


This is consistent with a qnark model interpretation of the composite nonet held: 




Qaq 


( 2 . 6 ) 


identifying qi,q 2 ,q 3 = u,d,s. Specihcally, Eq. ( p.6|) states that iV| is composed of one 
strange qnark and one strange antiqnark, k of one non-strange qnark and one strange an- 
tiqnark while Oq and {Nl + iV|)/v^ have zero strange content. Thns the ordering in Eq. 
( p.5|) natnrally follows if the strange qnark is heavier than the non-strange qnark, as has 
been well established. This ideal mixing pictnre works well for the vector mesons (with the 
rei dent ideations > 0, {Nl -|- A^'D/v^ —>• a;, K ^ K* and Oq —>■ p) and reasonably well 

for most of the other observed meson mnltiplets (see page 98 of [^). The exceptions are 
the low-lying 0“ and O’*" nonets. It is generally accepted that the deviation of the 0“ nonet 
from this pictnre can be nnderstood from the special connection of the psendoscalar havor 
singlet with the U{1)a anomaly of QCD. The case of the 0^ nonet has been less clear, in 
part becanse the existence of the scalar states needed to hll np a low-lying nonet has been 
difhcnlt to establish. 

Now a long time ago, Jaffe snggested that the low-lying scalars might have a qnark 
snbstrnctnre of the form qqqq rather than qq. This model can be pnt in the identical form as 


onr previons disenssion of Eqs. (|2.1|) - (|2.4| ) by introdneing the “dnal” havor qnarks (actnally 
diqnarks): 


Ta = eaicfif, T‘ = 


(2.7) 


wherein it shonld be noted that the qnark fields are anticommnting qnantities. Then we 
shonld write the scalar nonet as 






sdds sdus sdud 
suds situs suud 
adds udus udud 


( 2 . 8 ) 


In the present qqqq case both sides of Eq. (^^) shonld be taken to be positive. The 
tentative identiheations /o(980) = (A^j^-|-iV|)/-\/2 and a = wonld then lead to an ordering 
opposite to that of Eq. 


(/o) = (oo) > (o') • 


(2.9) 


This is in evident good agreement with the experimentally observed equality of the /o(980) 
and ao(980) masses. Furthermore it is seen that the ordering in Eq. ( p.9|) agrees with 
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the number of underlying (true) strange objects present in each meson according to the 
alternative ansatz ( p.8|) . 

If additional terms Q are added to the ideal mixing model in Eq. (O to yield 

^^mass = -ai:i{NN) - UIi{NNM) - cTi{N)Tr{N) - dTr(Ar)Tr(iVA^), (2.10) 


the states {Nl + A'"!) /a/ 2 and iV| will no longer be diagonal. The physical states will be 
some linear combination of these. This “non-ideally mixed” situation will be seen to be 
required in order to explain the experimental pattern of scalar decay modes. We would like 
to stress that, in the effective Lagrangian approach, no more than the assumption of mass 
terms like (|2.1(]|) is required; it is not necessary to assume a particular quark substructure 

), a linear combination 


for N^. That field may represent a structure like (|2.6|) , one like 
of these or something more complicated. Of course, it is still interesting to ask whether the 


resulting predictions are closer to those resulting from (^) or from 

A natural question concerns the plausibility of the “dual” ansatz in Eq. (p.8|) , which 
at hrst sight seems merely contrived to yield the ordering in Eq. (|2.9|) . Jaffe ||2^ showed 
that there is a dynamical basis for such an ansatz in the MIT bag model |^. It essentially 
arises from the strong binding energy in such a conhguration due to a hyperhne interaction 
Hamiltonian of the form 


■ SjF, ■ F 


( 2 . 11 ) 


where A is a positive quantity depending on the quark or antiquark wave functions. ^ ~ ^ 

is the spin operator and F = — (A are the Gell-Mann matrices) is the color-spin operator. 
The sum is to be taken over each pair (i, j) of objects (i.e. qq, qq or qq) in the hadron of 
interest. Eq. CT ) represents an approximation to the hyperhne interaction obtained from 


one gluon exchange in QCD; it is widely used in both quark model ||28| as well as bag model 
treatments of hadron spectroscopy. 

Standard application of ( p.ll|) to the p — n and A — N mass differences in the simple 
quark model yields: 


(tt] Tf/jj I tt) 

{N\H„\N) = -iA„„ 


iflHh/lp) = + 5 ^,,-, 
{A\H„\A} = +iA„„ 


( 2 . 12 ) 


*We are neglecting a possible term —eTr(AAf)Tr(AAf) which is second order in symmetry 
breaking. 
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in which a subscript has been given to the A factor for each quark configuration. It can 
be seen that A is expected to be fairly substantial - of order of several hundred MeV - in 
these cases. The evaluation of the expectation value of Eq. ( |2.11|) for the lowest scalar qqqq 
nonet state H is more complicated than for the above cases and yields a large enhancement 
factor due to the color and spin Clebsch-Gordon manipulations: 


{0+|ffv|0+) « -2.71A„55. 


(2.13) 


Thus, quark model arguments make plausible a strongly bound qqqq configuration. It 
should be remarked that the lowest lying O’*" nonet state in the quark model which diago¬ 


nalizes Eq. (|2.11|) is a particular linear combination of state 1 in which the qq pair is in a 
3 of color and is a spin singlet and state 2 in which the qq pair is in a color 6 and is a spin 
triplet: 


|0+) ^ 0.585|1) + 0.811|2). 


(2.14) 


A derivation of Eqs. (|2.13|) and (p.l4|) is given in Appendix A 


III. SCALAR NONET MIXINGS AND TRILINEAR COUPLINGS 


First let us consider the consequences of the generalized mass terms (|2.1U| ), which allow 
for arbitrary deviations from ideal mixing. The squared masses of the Oq and k, are read off 
as 


(oo) = 2a + 2b 
(k) = 2a -|- (1 -I- x) b. 


(3.1) 


Using the basis (^A|, ) > the mass squared matrix of the two iso-scalar mesons is 

also read off as 


2m^ (k) — (ao) -|- 2c -|- 2dx \f2 [2c -|- (1 -|- x) d] 


V2 [2c + (1 + x) d\ 


(ao) -I- 4c -I- dd 


(3.2) 


In obtaining this result Eqs. (|3.1|) were used to eliminate the parameters a and b. The 
physical isoscalar states and squared masses are to be obtained by diagonalizing this ma¬ 
trix. Notice that the four parameters a, b, c and d may be essentially traded for the four 
masses. We will take [|^ the strange to non-strange quark mass ratio x to be 20.5 for 
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definiteness. Then, up to a discrete ambiguity, the mixing angle between the two isoscalars 
will be predicted. 

It seems worthwhile to point out that the structure of our mass formulas provides con¬ 
straints on the allowed masses. To see this, note that the diagonalization of (|3.2|) yields the 
following quadratic equation for d = {1 — x) d: 


Qd^ — 8 (oo) — {k) 

+ 3m^ (a) (/o) — 6m^ (k) (oo) + (oo) — 5 {Arn^ (k) — (ao)) = 0, (3.3) 

where d = m? (a) + im? (/q) — 2in? {k) and we have eliminated c according to 6c = <5 — 
(4 + 2x) d. Here a and /o stand respectively for the lighter and heavier isoscalar particles. 
In order for d to be purely real, required at the present level of analysis, we must have 


(flo) “ (k) + 3m^ (ag) (cr) + (/q) + 9m^ (a) m? (/q) 

< 12m^ (oo) nn? (a) + (fo) . (3.4) 


Taking m (/q) = 980MeV and m (uq) = 983.5 MeV, according to [^, and m (a) = 550 MeV 
from PI we hnd that (|3.4|) limits the allowed range of m (k) to 


685 MeV < m (k) < 980 MeV. 


(3.5) 


It is encouraging that our recent study of ttK scattering (see also j^) yielded a value 
for m (k) of about 900MeV, within this range. 

The physical particles a and fo which diagonalize (|3.2|) are related to the basis states iV| 
and (Nl + Nf)/V2 by 


(T \ / cos6*s —sin6*s \ \ 

fo) \ sinOs cosOs ) \ ) 


(3.6) 


which dehnes the scalar mixing angle Og. Since Eq. (O) for d is quadratic we expect two 
different solutions for the pair (c, d) and hence for Og when we £x the four scalar masses 
m{ao), m{K), m{a) and m(/o). A numerical diagonalization for the choice m{K) ~ 900 
MeV as above yields the two possible solutions 


(а) 9g ^ -21° (3.7) 

(б) Og^ -89°. 

Solution (a) corresponds to a a particle which is mostly (presumably qqqq type) while 
solution (6) corresponds to a which is {Nl + V|)/v^ (i.e. qq type). We see that when 
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m, (GeV) 


FIG. 1. Scalar mixing angle solutions as functions of 


deviations from ideal mixing are allowed, the pattern of low lying scalar masses is by itself 
not sufficient to determine the quark substructure of the scalars. This statement is based 
on (p.lO|) which contains all terms at most linear in the mass spurion A4. 

For the complete allowed range of in Eq. ( |3.5| ) the two (“small” and “large”) mixing 
angle solutions are displayed in Fig. [^. Notice that the small angle solution is zero for 
itIk ~ 800 MeV; this is approximately where c = d = 0, which would correspond to the dual 
ideal mixing situation. In our convention — | < 6^8 < ^. 

Next let us consider the trilinear scalar-pseudoscalar-pseudoscalar interaction which is 
related to the main decay modes of the light scalar nonet states. We denote the matrix of 
pseudoscalar nonet helds by (x). The general Slf (3) flavor invariant Nc/x/) interaction is 
written as 


+ ^Tr (iV) Tr (d^c/>d^c/>) + GTr (Nd^c/>) Tr 
+ DTr (TV) Tr (d^c/>) Tr (d^c/>), (3.8) 

where A, B, C, D are four real constants. The derivatives of the pseudoscalars were intro¬ 
duced in order that (|3.8|) properly follows from a chiral invariant Lagrangian in which the 
held 0^ transforms non-linearly under axial transformations. The chiral aspect of our model 
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is largely irrelevant to the discussion in the present paper but, for completeness, will be 
briefly treated in Appendix B. 

Notice that the hrst term of (|3.8| ) may be rewritten as 

2ATr - ATr (N) Tr - 2ATr (iV0^0) Tr (0^0) 

+ ATr (N) Tr (0^0) Tr (0^0). (3.9) 


Thus, if desired, the complicated looking hrst term of ( |3.8| ) may be eliminated in favor of 
the most standard form Tr {Nd^(f)d^(f)). Our motivation for presenting it in the way shown is 
that, by itself, the hrst term of (|3.8|) predicts zero coupling constants for both /o —>■ tttt and 
a —^ KK when the “dual” ideal mixing identihcations, a = A| and /o = {Nl + A"|)/a/ 2, 
are made. This is in agreement with Jahe’s picture (see Section VB of ||2^) of the dominant 
scalar decays arising as the “falling apart” or “quark rearrangement” of their constituents. 
It is easy to see from that cannot fall apart into KK and that {Nl + A"|)/a/ 2 
cannot fall apart into tttt. 

Of course /o ^ tttt must be non-zero because /o(980) is observed in tttt scattering. In 
fact it also vanishes with just the term Tr (A^0^00^0) and the “conventional” identihcation 
a = {NI + N 2 )/\/2 and /o = N^. Our model contains two sources for /o — tttt: the deviation 
from ideal mixing due to the c and d terms in (|2.1CI|) and also the presence of more than one 
term in ( p.8|) . Note again that the use of ( p.lOj ) and (|3.8|) does not require us to make any 
commitment as to the quark substructure of N^. 

Using isotopic spin invariance, the trilinear Ncjxj) interaction resulting from (|3.8|) must 
have the form 


— d^N4)(t> — 


'^kK'k 


{pi^KT ■ df^TZK h.c.) 


+ l^ad,Kd,K + ■ d,7r + lI^f,d,Kd,K 


+ 


laoKK 


d^Kr ■ aod^K {i^df^Kd^rj + h.c.) + {i^d^Kd^rj' h.c.) 


+ 7ao7r77ao ■ df,7rdf,r] + yaoTr^ao ■ d^nd^rj' 

+ 'y<TrtriCrd^r]d^r] 'yan'r^'crd^T]’ 

+ Ifowfod/^vd/^V + lfovv'fodi,vdi^v’ + Ihv’v’fodi^v’d,,v', 


(3.10) 


where the y’s are the coupling constants. The helds which appear in this expression are the 
isomultiplets: 


K 



K = (^K- K^) , K = 
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(3.11) 


^ ^ (tTi T 17 ^ 2 ) , 

0 ~ (®oi T *002)) 


TT = TTa, 


Oq — ao3, 


in addition to the isosinglets a, /o, r] and rj'. The expressions for the 7 ’s in terms of the 
parameters A, B, C and D as well as the scalar and pseudoscalar mixing angles are listed, 
together with some related material, in Appendix C. Notice that if we restrict attention to 


those terms in which neither an rj nor an rj' appear [hrst six terms of (|3.10|) ], their coupling 
constants only involve two parameters A and B. These are the terms which will be needed 
for the subsequent work in the present paper. 


IV. TESTING THE MODEL’S COUPLING CONSTANT PREDICTIONS 

Now let us consider how well the hve coupling constants '')k,k-ki la-wir, laKK, Ifomr and 
IfoKK, can be correlated in terms of the two parameters A and B. These coupling constants, 
which are listed in Eqs. OS are the ones which are relevant for the discussions of tttt 
scattering given in and tiK scattering given in [^]. 

A very important question concerns the way in which these 7 ’s are to be related to 
experiment. For an “isolated” narrow resonance the magnitude of the coupling constant is 
proportional to the square root of the width. Actually, the only one of the hve for which 
this prescription roughly applies is ’jfoTrn', the appropriate formula is given in Eq. (4.5) of [0. 
Even here there is a practical ambiguity in that, while the tttt branching ratio is listed in 
, the total width is uncertain in the range 40 —lOOMeV. The determination | 7 /o 7 r 7 r| = 2.43 
GeV“^ given in p| is based on using Ttot (/o) as a parameter in the model analysis of tttt 
scattering and making a best ht. 

The situation for ’J/oKK is somewhat similar due to the poorly determined rtoi(/o)- 
There is an additional difficulty since the central value of the /o(980) mass is below the KK 
threshold. Thus the value | 7 /Q_ft:x| ~ 10 GeV“^ presented in Section V of p, is based on a 
model taking the hnite width of the initial state into account. Incidentally, the non-negligible 
branching ratio for /o KK in spite of the unfavorable phase space is an indication that 
the /o “wavefunction” has an important piece containing ss. 

The a, as “seen” from the analysis of @], for example, is neither isolated nor narrow. A 
suitable regularization of the tree amplitude near the a pole was argued to be of the form: 

m^G m„G 

— s — s — 
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coupling constant 

value (GeV 

obtained from 

IT/ottvtI 

2.4 

TTTr scattering 

hfoKKl 

« 10 

TTTr scattering 

iTcrTTTT 1 

7.8 

TTTT scattering 

\'yK.K7r\ 

5.0 

ttK scattering 

luKK 

8 

■kK scattering 


TABLE I. Coupling constants previously obtained in [2] and [3]. 


The situation concerning 'y^K-w is similar to the one for 
to the / = - amplitude of ttK scattering (see Section IV of 0) yields lyK/^Trl 

values obtained as 


where G and G' are real. G is taken to be proportional to while G' is considered to 
be a regularization parameter. For a narrow resonance with negligible background it would 
be expected that G' = G. However, considering both G and G' as quantities to be ht (or 
essentially equivalently, restoring local unitarity in a crossing symmetric way) yields G' ^ G. 
The determination | 7 o- 7 r 7 r| = 7.81 GeV“^ is based on such a ht. 

Tjr. Making an analogous ht 
5 GeV^ 

This value, however, is based on inputting the | 7 /o 7 r 7 r|, \lfoKK\ and ly^ 
above and making a particular choice of jaKK- The value of 'jaKK was however not very 
accurately determined in this model; a compromise choice was 'J^kk ~ 8 GeV“^. 

A summary of the coupling constants previously obtained is shown in Table |. 

The discussion above illustrates that it seems necessary to obtain the coupling constants 
of the low-lying scalars from a detailed consideration of the relevant scattering processes. It is 
not sufficient to read them oh from at the present time. Furthermore their interpretation 
is linked to the dynamical model from which they are obtained. 

It seems to us that a relatively clean way to test the correlation between the coupling 
constants in Table is to recalculate the nK scattering amplitude and, instead of taking 
\lfo-rrn\, \lfoKK\ and lyo-TTTrl from the tttt scattering output and regarding ■jKKn and ■yaKK as 
htting parameters as in |^, just A and B are now taken to be htted. 

We work within the same theoretical framework that was developed in for the tttt 
scattering analysis and was further explored in for the case of nK scattering. In this 
framework, the ttK scattering amplitude is computed in a model motivated by the l/W 
picture of QGD and its real part is given as a sum of regularized tree level graphs which 
include all resonances that contribute to the amplitude up to the energy region of interest. 
The relevant Feynman diagrams are shown in Fig. 1 of 0. 

In the I = - channel, we perform a fit) using the MINUIT package, of this model to 
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the experimental data. Specifically, in addition to A and B, the parameters to be fit are the 
regularization parameter in the k propagator, G'^ (which can also be interpreted as a total k 
decay width), and parameters of the resonance it'Q(1430): its mass M*, its coupling 7 * and 
the regularization parameter in its s-channel propagator This will be done for different 
choices of Note that the scalar mixing angle Og (see Section III) will be different for each 
choice of In fact, as already discussed, this actually gives two different mixing angles for 
each m^, one (large angle solution) closer to the qq ansatz (| 2 . 6 |) and the other (small angle 
solution) closer to the qqqq ansatz (|2.8|) . It is very interesting to see which one is chosen in 
our model. More details of the model are given in [Q. The possible values of are limited 
by (p.5|) for consistency with our present model for masses based on Eq. (|2.10|) . 

Let us first choose = 897MeV, as obtained in |^. Then the fit [] to the real part 
1 1 

of the / = - amplitude, Rg is shown in Fig. while the fitted parameters and resulting 

predicted coupling constants are given in Table |I|. The results for both possible mixing 

1 

angles corresponding to = 897MeV are included. It is seen that the fits to Rg are 
essentially equally good compared to each other and compared to the one in 0 . However if 
we compare the coupling constants in Table with those obtained previously in Table | we 
see that while the coupling constants XfoTr-K, 'IfoKK, lann and 'JkKtt obtained with 9s ~ — 20 ° 
agree with those obtained earlier in connection with tttt and nK scattering, their values 
obtained with 9s ~ —89° do not agree so well. 

Furthermore the value of 'yfomr obtained with 9s ~ —89° would lead to a value for the 
/g width several times larger than the experimentally allowed range. It thus seems that the 
qqqq picture, to which 9s ~ —17° is much closer, gives a better overall description of the 
scalar nonet than does the qq picture. 

It is interesting to investigate the effect of changing within the range given in Eq. 
(fj.5|). As examples. Tables and show the fitted parameters for = 875 MeV and 

= 800 MeV respectively. Several trends can be discerned. As decreases from 897 
MeV the goodness of fit actually improves from = 3.94 to = 2.3 at = 800 MeV. On 
the other hand the value of y/oTr-n- increases so that at = 875 MeV the /g —>• tttt width is 
in slightly better agreement with experiment and at = 800 MeV it is many times larger 
than allowed by experiment. It seems that the fit at = 875 MeV is not very different 
from the one at = 897 MeV; this gives an estimate of the “theoretical uncertainty” in 
our calculation. On the other hand = 800 MeV seems to be ruled out, as are still lower 


iThe experimental data points are taken from ||30t|. 
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values of 

3 

Another argument in favor of the larger values of m,^ can be made by examining the / = - 
7tK amplitude shown in Fig. It is seen that decreasing worsens the agreement with 
experiment. This feature arises because jaKK, to which the / = - amplitude is sensitive, 
increases with decreasing This situation was discussed in more detail in section V of 
where it was noted that higher mass resonances may be important in this channel. 

We note that the three parameters describing the A’g(1430) are stable to varying m^. 
All the hts yield for the parameters A and B that — 1. Using (|3.8|) then shows that 
approximately looks like 


^ 2A [Tr (Nd^cpd^cP) - pTr (N) Tr (9^09^0)] + 


(4.2) 


where p is a positive number slightly less than unity and the three dots stand for the C and 
D terms which only contribute to vertices involving at least one p or p'. 

Using this model we can also estimate the partial decay width of ao(980) —^ KK which 
is entirely determined in terms of the parameter A [see Eq. (|C4|) ]. As in the case of /o(980), 
the resonance lies below the decay threshold so the effect of the hnite width of the decaying 
state must be taken into account [see for example footnote 2 of [0|]. The results are shown 
in Table V (taking = 897 MeV) corresponding to the extremes of the total width range 
given in Also the effect of the mass difference between the charged and neutral kaons 
is taken into account. The numerical values seem reasonable. 


V. DISCUSSION 


We studied the family relationship of a possible scalar nonet composed of the /o(980), 
the ao(980) and the a and k type states found in recent treatments of mi scattering and 
tiK scattering. The investigation was carried out in the effective Lagrangian framework, 
starting from the notion of “ideal mixing”. First it was observed that Okubo’s original 
treatment allows two solutions: one the conventional (e.g. vector meson) qq type and the 
other a “dual” picture which is equivalent to Jaffe’s qqqq model. 

The four masses of our scalar nonet candidates have a similar, but not identical pattern 
to the one expected in the dual ideal mixing picture. In order to allow for a deviation from 
ideal mixing, we have added more terms to the Lagrangian [see ([^.1U|)]. The resulting mass. 


iThe experimental data points are taken from |31|. 
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Vs (GeV) 

1 /2 

FIG. 2. Comparison of the theoretical prediction of Rq with its experimental data. 



3/2 

FIG. 3. Comparison of the theoretical predictions of Rq with its experimental data. 
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Eitted Parameter 

Os = -20.33 

Os = -89.14 


314 ± 3 MeV 

322 ± 3 MeV 


1390 ± 4 MeV 

1389 ± 4 MeV 

7 * 

4.42 ± 0.09 GeV^ 

4.4 ±0.09 GeV^ 


275 ± 10 MeV 

273 ± 11 MeV 

A 

2.51 ±0.03 GeV^ 

2.57 ±0.03 GeV^ 

B 

-1.95 ±0.04 GeV^ 

-2.12 ±0.04 GeV^ 


3.94 

3.95 

Predicted Couplings 

'^kK'K 

-5.02 GeV^ 

-5.14 GeV^ 

O^CTTTTT 

7.26 GeV^ 

4.33 GeV^ 

T/tttt 

1.46 GeV^ 

-6.56 GeV^ 

laKK 

9.62 GeV^ 

13.69 GeV^ 

IfKK 

10.10 GeV^ 

-5.78 GeV^ 


TABLE II. Extracted parameters from a fit to the nK data, = 897 MeV. 


Eitted Parameter 

= -15.61 

m 

OO 

C7) 

1-^ 


346 ± 2 MeV 

357 ± 3 MeV 

M, 

1389 ± 4 MeV 

1388 ± 4 MeV 

7* 

4.42 ± 0.09 GeV^ 

4.39 ± 0.09 GeV^ 

g: 

275 ± 10 MeV 

272 ± 10 MeV 

A 

2.87 ±0.03 GeV^ 

2.96 ± 0.03 GeV^ 

B 

-2.34 ±0.03 GeV^ 

-2.56 ±0.04 GeV^ 


3.23 

3.26 

Predicted Gouplings 

'y kKtv 

-5.75 GeV^ 

-5.92 GeV^ 

O^CTTTTr 

8.36 GeV^ 

-4.58 GeV^ 

T/tttt 

2.53 GeV^ 

8.13 GeV^ 

laKK 

10.45 GeV^ 

-15.62 GeV^ 

IfKK 

12.76 GeV^ 

8.30 GeV^ 


TABLE III. Extracted parameters from a fit to the nK data, = 875 MeV. 
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Fitted Parameter 

= -0.84 

0s = 71.37 


450 ± 2 MeV 

479 ± 2 MeV 


1387 ± 4 MeV 

1384 ± 4 MeV 

7* 

4.40 ± 0.09 GeV^ 

4.36 ± 0.09 GeV^ 

g: 

273 ± 10 MeV 

268 ± 11 MeV 

A 

4.32 ± 0.03 GeV^ 

4.50 ± 0.04 GeV^ 

B 

-3.91 ±0.03 GeV^ 

-4.29 ± 0.04 GeV^ 


2.34 

2.39 

Predicted Couplings 

'ynKiv 

-8.64 GeV^ 

-9.01 GeV^ 

'^U'K'K 

11.76 GeV^ 

-4.15 GeV^ 

T/'TTTT 

7.65 GeV^ 

14.52 GeV^ 

luKK 

11.41 GeV^ 

-20.91 GeV^ 

IfKK 

24.12 GeV^ 

19.85 GeV^ 


TABLE IV. Extracted parameters from a fit to the nK data, = 800 MeV. 


decay widths 

r“ = 50 MeV 

r“ = 100 MeV 

r(aO- 

o 

t 

0.924 MeV 

2.049 MeV 

1 

oo 

l-H 

K+K-) 

1.371 MeV 

2.455 MeV 


TABLE V. Predicted Qq ^ KK decay widths 
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mixing and scalar-pseudoscalar-pseudoscalar coupling patterns [see (|3.8|) ] were discussed 
in detail. The outcome of this analysis is that the dual picture is in fact favored. More 
quantitatively, the appropriate scalar mixing angle in Eq. (|3.6| ) comes out to be about 
—17° ±4° compared with 0° for dual ideal mixing and ±90° for conventional ideal mixing. 
This corresponds to ranging from 865 — 900 MeV. 

The coupling constant results obtained here may be useful for a number of applications in 
low energy hadron phenomenology. These are dehned in Eq. ( |3.10|) and listed in Appendix 
C. Typical values of A and B may be read from the small magnitude angle solution in Tables 
iand|Tl|. We expect to improve and further check the accuracy of this model by extending 
the underlying models of tttt and ttK scattering to higher energies and to other channels. 
Finally, it may be interesting to compare our results with those of quark model and lattice 
gauge theory approaches to QCD. 


ACKNOWLEDGMENTS 

The work of D.B., A.H.F. and J.S. has been supported in part by the US DOE under 
contract DE-FG-02-85ER 40231. The work of F.S. has been partially supported by the US 
DOE under contract DE-FG-02-92ER-40704. 


APPENDIX A: DIAGONALIZATION OF HYPERFINE HAMILTONIAN 


In this Appendix, we give some explicit details of the derivation of (p2.13| ) and (p2.14| ) 
which, while not being explicitly used in our approach, furnish the main reason for expecting 
the scalar qqqq states to be especially strongly bound. Our results agree with those of Jaffe 
who followed a different method. 

Let us begin by considering only flavor quantum numbers in order to write down the 
quark content of members of a qqqq scalar nonet. Taking the quarks to be in the fundamental 
representation, 3, of SU (3)/ we have the familiar irreducible decomposition of products of 
quark states: 


3(g)3 = 603 (Al) 

3(8)3 = 603. (A2) 

So the only possibility for obtaining a qqqq flavor nonet is from the combination 3 3 of 

(g) states. Let g* be a basis for the representation space 3, where i=l,2 and 3 correspond 
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to up, down and strange quarks respectively, with conjugate (antiquark) basis (p. Then we 
can consider “dual quark” bases corresponding to the qq and qq flavor triplet spaces (thus the 
states are antisymmetric with respect to exchange of flavor indices), namely := emjkQ^q^ 
and T™ := e"^^^qjqk- Up to (anti)symmeterization and linear combinations we have the 
flavor nonet given in Eq. (|2.8|) . Since and contain at most one strange quark each 
the nonet states contain at most two strange quarks. We note also that, in contrast to the 
conventional qq scalar nonet, is non-strange in this realization. 

In order to complete the description of qqqq scalar nonets we consider the spin and color 
quantum numbers. Using the facts that (i) the qq and qq parts of the state are individually 
totally antisymmetric and (ii) the overall qqqq hadron must be a color singlet, where the 
quarks transform according to the fundamental representation of SU{3)c, we obtain just 
two possibilities which include scalar flavor nonets (noting that 6c (8) 6c = Ic 08c 0 27c), 
namely 


| 0 +, 9 ), 


0+,3f,3, 


0+,3f,3c 

qq L 


(A3) 


lot 9)2 


U,3,,6< 


J qq 


l+,3f,6. 


_? 

qq 


(A4) 


where we have shown the spin-parity, flavor and color representations respectively for qq and 
qq separately. 

The “hyperflne” interaction Hamiltonian needed for our discussion is given in Eq. ( p.ll ). 
Given two representations of SU(n) we have the well-known relationship between the 
quadratic Casimirs of these representations, say and and that of their product: 


J A ■ 3 B — — 


2 L 


* total 


'B 


(AS) 


It can be seen, using (^^), that the parts of the hyperflne Hamiltonian which involve 
sums over qq or qq pairs are diagonal with respect to the bases for the scalar nonets chosen 
in (^) and (|A4|) . In order to calculate the expectation value of the qq terms in (p.ll|) using 
( |X5|) we first expand the bases (|A3|) and (|A4|) in terms of states where the spin and color of 
the qq pairs are explicit. 

To And the recoupling coefficients we follow Close |^|, where more detail is given. For 
the case of spin recoupling we have, assuming that all of the quarks in the scalar meson are 
in relative s-wave states, that in order to couple to total angular momentum J = 0, either 
both qq pairs must be in = 1~ or both in = 0“ states, which we denote as vector, 
(U), and pseudoscalar, (P) respectively. Thus we can expand the spin part of the state in 
the following manner: 
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\Jtotal — 0)lor2 ~ OiPP + (3VV, 


(A6) 


where a and (3 can be determined in each case by rewriting both sides (the left-hand-side 
will be different for the two states (|A3|) and ([A4| )) in terms of their constitnent qnarks and 
antiqnarks nsing the usnal Clebsch-Gordon identities for SU{2). 

Similarly for the color states we note that, since 3(^3 = 8 01, only combinations of 
the form 


|8e)q,-(g) \8c)qq + (3’\lc)qq<S> m (^7) 


include color singlets and therefore the color parts of (|A3|) and (|A4|) can be written in 
terms of this basis. For brevity we simply present the results of our recoupling coefficient 
expansions in Table 


nonet 

spins oiqq pairs 

color products oiqq pairs 

|0+,9)i 

|0+,9)2 

ipp + ^VV 
^pp - Iff 

■^lc{8)lc - Y^8c{g)8c 
Y^lc (g) Ic + y^8c 0 8c 

TABLE VI. Spin and color recouplings for flavour nonets 


Representation 

p 2 

3 or 3 

4 

3 

8 

3 

1 

0 

6 

10 

3 


TABLE VII. SU(3) Quadratic Casimirs 


In order to give an idea of the next step let us look at one of the off-diagonal elements of 
{Hfif), where H^f is as in ( |2.11| ), with respect to the basis given in (|A3|) and (|A^). Labelling 
the quarks/antiquarks qiq 2 Q 3 Q 4 : we have that the only non-vanishing off-diagonal pieces in 
(Hhf) are the sums over (13), (14), (23) and (24). For example, applying (|A5|) yields 


Si-S3Fi-F3|0+,9f)i = 


-Mpp + l.^vv 

4 2 2 2 


3 ' v/3 


1 . - 


: 8 c 


I 8c 


(A8) 
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where for the color operators we have used the SU{3) Casimirs given in Table [VII| . Finally 
we take the inner product with the expansion of 9)2 in Table |V| which gives that 


(Si ■ S 3 F 1 ■ F3)21 


1 

4 



(A9) 


There are, as noted above, four such combinations, all of which contribute equally. An 
analogous calculation can be performed for the diagonal matrix elements giving hnally: 


{Hhf)ab — —A 




(AlO) 


where a and b run over the indices 1 and 2 labelling the flavor nonets |0^, 9)^ and |0^, 9 ) 2 - 
Thus the eigenstates of the hyperhne interaction correspond to mixtures of these nonets, 
corresponding to energies Ei = —2.71A and E 2 = —0.12A, which are in agreement with 
m- The corresponding eigenstates are: 


|0+, 9) = 0.585|0+, 9)1 + 0.811|0+, 9)2 

|0+, 9*) = 0.811|0+, 9)1 - 0.585|0+, 9 ) 2 - (All) 


APPENDIX B: CHIRAL COVARIANT FORM 

Here we present the terms of the total Lagrangian involving the scalar nonet A'"^(x) in 
chiral invariant or (for the mass terms which break the chiral symmetry) in chiral covariant 
form. We follow the general method of non-linear realization described in but our 
notation is as in Appendix B of P]. The object ^ = exp{i(j)/E t,) discussed there transforms 
as 


e ^ = KiUi (Bl) 

under chiral transformation. Our nonet held is considered to transform as if it were made 
of “constituent” quarks, namely 

N KNKl (B2) 

With the convenient objects 

we write the additional Lagrangian terms involving N: 
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C = -^Ti iV^NV^N) - aTr(iViV) - ^Tr 


'nN - cTr(A)Tr(A) 


^Tr(A)Tr [n 
CTt{Np^)Tt{p^) + DTt{N)Tt{p^)Tt{p^)] 


Ae^^^e,efNi{p,)l{p,)i + 5Tr(A)Tr(p^p^) 

(B4) 


where V = — iv^ and M. = AlMs the spur ion matrix dehned after (|2.1| ). The entire 

Eg■ (p^) is formally invariant if we allow Af to transform as A1 —> f/iAlf/jj. This La- 
grangian reproduces (|2.10|) and (|3.8|) but also contain interactions with extra pions. These 
extra interactions do not change anything in this paper or in the tree-level formulas for (jycf) 
scattering in and |^. 


APPENDIX C: COUPLING CONSTANTS 


Here we hnd the scalar-pseudoscalar-pseudoscalar coupling constants dehned in (|3.10|) 
in terms of the parameters A, B, U, D [see ( |3.^ )], the scalar mixing angle [see (pT6|)] and the 
pseudoscalar mixing angle, Op. The latter is dehned according to: 



cos^p —sin^p 
sin^p cosOp 


(</>l + </> 2 )/V 2 


(Cl) 


where p and p' are the helds which diagonalize the pseudoscalar analog of (^). The usual 
convention employs a diherent basis; in this convention the angle is 6^ and 



cosOu -sin6'„\ / {(j)\ + (j)l - 2(j)l)/^/Q\ 

sinOu cosOu ) \ {(j)\ + (j)l + (j)l)/VQ ) 


(C2) 


The relation between the two angles is 


9p = 9u + 54.74° ^ 37° (C3) 

in which case (see for example 9u ~ —18° was taken. More recent analyses ( and 
1 ^ ) have modihed this treatment somewhat by considering derivative mixing terms as well 
as non-derivative ones. 

Note that the basis for ( |C1|) was chosen so that qq is the more natural picture for the 
pseudoscalar nonet, in contrast to (|3.6|) for the scalars. Because the mixing angles can take 
on any values, this in no way biases the analysis one way or the other. 

The y’s are predicted in the present model as 
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IkKtt = 'laoKK = —2A 

7 ^^^ = 2Bsm9s — V2{B — A)cos9s 
'faKK = 2(25 — yl)sin6's — 2\/2Bcos9s 
T/otttt = V2{A — B)sm9s — 2Bcos9s 
IfoKK = 2(A — 2B)cos9s — 2\/2Bsm9s 
7 kKt] = Csm9p — V2{C — A)cos9p 
= y/2(A — C)sin9p — Ccos9p 
lao-KT) = {C — 2A)sm.9p — \/2Ccos9p 
lao-Ki)' = (2A — C)cos9p — \/2Csm.9p 

= V2{B + D) — ^{C + 2A + 4D)sin29p + V2{C + D)cos^9p sin9s 

— {B + D) - + 2D)sm29p + {A + D)cos^9p + Csm'^9p cos9s 

v2 

lari'ri'= V^(5 + 5) + -(C + 2A + 45)sin26*p + v^(C + 5)sin^6'p sm9s 

— {B + D) + + 2D)sm29p + {A + 0)8111^9p + Ceos‘S9p cos9s 

v2 

7 ( 7»?7 = [v^(C' + 5)sin26'p + (C + 2^4 + 45)cos26'p sin 6*5 

— |^V^(C + 2D)cos29p + [A — C + 5)sin26'pj co89s 

'yfonp= —V2{B + D) + ^{C + 2A + 4:D)sm29p — V2{C + D)cos^9p cos9s 

— (5 + 5)- + 2D)sm29p + {A + D)cos^9p + Csm'^9p sm9s 

v2 

= — v^(5 + 5) +-(C + 2^4 + 45)sin26*p + v^(C + 5)sin^6'p cos9s 

— (5 + 5) + ^=(C + 25)sin26*p + (A + 5)sin^6*p + Ccos^ 6 *p sm9s 

v2 

7/oW = “ V2{C + D)sm29p + {C + 2A + 4D)co829p co89s 

— \f2{C + 2D)cos29p + [A — C + 5)sin26'p sm9s 


(C4) 

(C5) 

(C6) 

(C7) 

(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 


(C14) 

(CIS) 


(C16) 


(C17) 

(C18) 
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